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Abstract 

(N 

, , We construct a two-dimensional diffusion process with rank-dependent local drift and disper- 

sion coefficients, and with a full range of patterns of behavior upon collision that range from 
Oh totally frictionless interaction, to elastic collision, to perfect reflection of one particle on the 

other. These interactions are governed by the left- and right-local times at the origin for the 
distance between the two particles. We realize this diffusion in terms of appropriate, appar- 
ently novel systems of stochastic differential equations involving local times, which we show 
are well posed. Questions of pathwise uniqueness and strength are also discussed for these 
^ systems. 

The analysis depends crucially on properties of a skew Brownian motion with two- valued 
drift of the bang-bang type, which we also study in some detail. These properties allow us to 
^ compute the transition probabilities of the original planar diffusion, and to study its behavior 

under time reversal. 

o 
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1 Introduction 



We construct a planar diffusion (Xi(-), X 2 (-)) according to the following recipe: each of its com- 
ponent particles Xl(-) and X 2 (-) behaves locally like Brownian motion. The characteristics of 
these random motions are assigned not by name, but by rank: the leader is assigned drift —h < 
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and dispersion p > , whereas the laggard is assigned drift g > and dispersion a > . One 
of the dispersions is allowed to vanish, but not both; similarly for the drifts. In the interest of 
concreteness and simplicity, we shall set 



A := g + h > , p z + a z = 1 . (1.1) 

A bit more precisely, we shall construct a complete probability space 5, P) endowed with a 
filtration F = {#(t)}o<t<oo that satisfies the "usual conditions" of right continuity and of augmen- 
tation by P— negligible sets, and on it two pairs (-E>i(-) ; B 2 (-)) and (Xi(-) , X 2 (-)) of continuous, 
F— adapted processes, such that (Bi(-), B 2 (-)) is planar Brownian motion and (Xi(-), X 2 (-)) a 
continuous planar semimartingale that starts at some given site (Xi(0), X 2 (0)) = (x\,x 2 ) G M 2 
on the plane and satisfies the dynamics 

dXi(t) = (g l{xi(t)<xi(t)} - h l{Xi(t)>x 2 (t)}J dt + (pl{Xi(t)>x 2 (t)} + <? l{Xi(t)<x 2 (t)}J dB x {t) 

+ dL x ^(t) + dL" 1 (t) , (1.2) 

dX 2 (i) = fpl^^^W} -^ 1 {x 1 (t)<x 2 (t)}') dt + rpl{Xi(t)<x 3 (i)} + °" l{x 1 (t)>x 2 (t)}J d5 2 (t) 

+ l^ d L Xl -^(t) + l^dL^-^(t) . (1.3) 

Here and in the sequel we denote by L x {-) = L x {- ; 0) the right-continuous local time accu- 
mulated at the origin by a generic continuous semimartingale X(-) , by L x {-) = L~ x (- ; 0) its 
left-continuous version, and by L x (-) = (L x (-) + L x (-))/2 its symmetric version; we collect in 
section 2 the necessary reminders from the theory of semimartingale local time. 

Each time the particles collide, their trajectories are "dragged" by amounts proportional to the 
right local times, L Xl ~ X2 {t) and L X2 ~ Xl {t) , respectively, that have been accumulated up to that 
instant t at the origin by the differences Xi(-) — X 2 (-) and X 2 (-) — Xi(-) ; this is the significance 
of the last two terms in each of (1.2), (1.3). With the notation 

C:=l + ^. V-=l-^, (1-4) 

the proportionality constants of these interactions, Q and r]i for Xj(-) (i = 1, 2), will be assumed 
to satisfy the conditions 

C + r^^O, <1. (1.5) 

v + C 

We shall discuss in detail the significance of these conditions for the system of equations 
(1.2), (1.3); in particular, the fact that they are not only sufficient but also necessary for the well- 
posedness of the above system of stochastic equations. For the time being, let us note that in the 
special case d = rji — 1 , the trajectory Xl(-) of the first particle crosses the trajectory X 2 (-) 
of the second particle without "feeling" it, that is, without being subjected to any local time drag; 
as we shall see in subsection 5.2, we obtain this same effect under the more general condition 
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(5.5). Likewise, the second particle does not "feel" the first, when ( 2 = r/ 2 = 1 or, more gener- 
ally, under the condition (5.6). When Q = r/i = 1 (i — 1, 2), the local times vanish completely 
from (1.2), (1.3) and we are in the situation studied in detail by Fernholz et al. Fernholz, 
E.R., Ichiba, T., Karatzas, I. & Prokaj, V. (2012) (2012). In this case, the collisions of the 
particles are totaly frictionless. 

At the other extreme ( = ^ r] (respectively, i] = ^ () the trajectory X\(-) of the first 
particle bounces off the trajectory X 2 {-) of the second particle (resp., the other way round), as if 
the latter trajectory were a perfectly reflecting boundary; cf. subsection 5.1. Think of the second 
(resp., the first) particle as being "heavy", so that in collisions with the "light" first (resp., second) 
particle its motion is unaffected, while the light particle undergoes perfect reflection. 

In between, for other values of the parameters, we have collisions that are neither totally fric- 
tionless (without local time drag), nor perfectly reflecting, but "elastic": The particles are subjected 
in general to local time drag, and this kind of friction manifests itself in an asymmetric fashion - 
due to the presence of both right- and left- local times at the origin L Y (■) = L Y (■ ; 0+) and 
L Y (-) = L Y (- ; 0-) in (1.2), (1.3) for the difference Y(-) = X x {) - X 2 {-) . We call such colli- 
sions "skew-elastic". 

1.1 Preview 

Under the conditions of (1.5), the system of equations (1.2), (1.3) will be shown in section 4 to 
admit a weak solution, which is unique in the sense of the probability distribution; cf. Theorem 
4.1. Using a common terminology: under the conditions of (1.5), the system of equations (1.2), 
(1.3) is well posed. 

A crucial role in establishing this result will be played by the properties of the difference 
Y(-) = Xi(-) — X 2 (-) , for which we show that 

W(t) := Y(t) - (a?i - x 2 ) + A f sgn(Y(s)) ds - 2 (2 a - l) L Y (t) , < t < oo 

Jo 

is standard Brownian motion W(-) in the notation of (1.1), (1.5). To put a little differently: we 
identify this difference Y(-) = Xi(-) — X 2 (-) as a so-called Skew Brownian Motion with Bang- 
Bang drift, a process that we study in detail in Section 6. 

Similarly, recalling the notation of (1.4) and setting 

a = ry (Ci + C2) + C(^7i + ^72) 

we identify the process 

V(t) := Xi{t) + X 2 (t) - (xi + X2 )-{g-h)t-2{l-P) L Y (t) , < t < 00 
as another standard Brownian motion, whose cross-variation with the Brownian motion W(-) is 

(V,W){-) = (V,Y)(.) = 7 / sgn(y(i))dt, 7 := p 2 -a 2 . 

Jo 
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These identifications allow us then to represent the motions X\(-) , X 2 (-) of the individual parti- 
cles in the form 

Xt{t) = x 1 +p l t + p 2 {Y+(t) - y+) - a 2 (Y-(t) - y~) + (l - - j) L Y (t) + paQ(t) , 

X 2 (t) = X2 + flt -a 2 (Y + (t) - y + ) + p 2 (Y-(t) - y~) + (l - /3 - 7) L Y (t) + paQ(t) ; 

here Q(-) is yet another standard Brownian motion, independent of the difference Y(-) — X±(-) — 
X 2 (-),and 

2 7 2 
p = g p — ha . 

This way we construct a weak solution to the system of equations (1.2), (1.3), and also show that 
uniqueness in distribution holds for it. 

Always under the conditions of (1.5), the system of equations (1.2), (1.3) is shown in section 4 
actually to admit a pathwise unique, strong solution; cf. Theorem 4.2. Here we refine the Le Gall 
LeGall, J.F. (1983) (1983) LeGall, J.F. (1984) (1984) methodology, that we used in the 
recent work Fernholz et al. Fernholz, E.R., Ichiba, T., Karatzas, I. & Prokaj, V. 
(2012) (2012) to establish pathwise uniqueness for a generalization of the perturbed TANAKA 
equation of Prokaj Prokaj (2010) (2011). 

In fact, the conditions in (1.5) turn out to be not just sufficient but also necessary for the well- 
posedness of the system (1.2), (1.3); cf. Proposition 6.1. As we shall see in Remarks 3.1 and 3.2, 
this system admits no solution in the case r] = —( 7^ ; whereas it has lots of solutions, i.e., 
uniqueness in distribution fails for the system of equations (1.2) and (1.3), when 77 = £ = . 
Finally, if we do have r] + ( ^ yet (1.5) fails because a ^ [0, 1] , it is seen in Remark 6.1 that 
the system (1.2), (1.3) once again fails to admit a solution. 

Section 5 discusses some special configurations of the parameters r] u Q (i = 1,2) in (1.2), 
(1.3) that give rise to some rather interesting structure. We see, in particular in the non-degenerate 
case pa > 0, that when (3 = (respectively, (3 = 2), the trajectory Xl(-) V X 2 (-) of the 
"leader" (respectively, Xi(-) A X 2 (-) of the "laggard") is Brownian motion with drift, with perfect 
reflection on the trajectory Xl(-) A X 2 (-) of the "laggard" (respectively, Xi(-) V X 2 (-) of the 
"leader"), which is then another, independent Brownian motion with drift. 

Section 6 develops the theory and properties of the Skew Brownian Motion with Bang-Bang 
drift. Finally, Section 7 uses these properties to compute the transition probabilities and the time- 
reversal of the planar diffusion (Xi(-), X 2 (-)) . 

1.2 Extant Work and Open Questions 

The study of multidimensional stochastic differential equations that involve a local time sup- 
ported on a smooth hypersurface starts with the work of ANULOVA ANULOVA (1980) (1980) 
and Portenko Portenko, N.I. (1976) (1976) Portenko, N.I. (1979) (1979) Portenko, 
N.I. (1979.b) (1979) Portenko, N.I. (1990) (1990). To the best of our knowledge, systems of 
stochastic equations of the type 

Xi(-) = x l + B l (-) + Y J QijL Xi - x '(-), i = l,---,n (1.6) 
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for a suitable array of real constants (%)i<ij< n , with Bi(-), - ■ ■ , B n {-) independent standard 
Brownian motions, were studied first by Sznitman & Varadhan Sznitman, A.S . & Varad- 
HAN, S .R.S . (1986). In fact, these authors consider the more general model 

N 

X(t) = x + B{t) + ^2q k L nk - x {t), 0<t<oo, (1.7) 
fc=i 

where X(-) : = (X x (•),•••, X n (-))' , B(-) := (-Bi(-), ■ ■ ■ ,B n (-))' is Brownian motion in M n , 
x G K n , the unit column vectors n k generate pairwise distinct hyperplanes, and the column 
vectors q k satisfy the orthogonality conditions q k ■ n k — for k — 1, . . . , N . When g = 
h = and a = p , it can be verified - using the relationships (3.16) between the symmetric 
local time and the right local time - that the system (1 .2)-(l .3) is equivalent to the model (1 .7) with 
parameters n = 2 , N = 1 , 

m := (1,-1)7 V2, and Ql .= (a(l -&) + (1 -^)(1 -a) , a{\ -( 2 ) + (1 -^)(1 -a))' . 

The orthogonality conditions amount then to = ( . Thus, we can apply the results of Sznitman 
& Varadhan Sznitman, A.S. & Varadhan, S.R.S. (1986), ifg = h = 0,a = p, 
i] = ( in our system (1.2)-(1.3). 

There are rather obvious similarities, as well as differences, between the system (1.6) and that 
of (1.2), (1.3). In particular, it would be very interesting to extend the results of this paper to 
systems of stochastic differential equations of the type 

n » . n „ . 

fc=i Jo k=i Jo 

+ E {4L Xi ~ Xi (-) + H j L Xj ~ Xi (-)) , z = 1, • • • ,n (1.8) 

for an arbitrary number n E N of particles, with xi, • • • ,x n and Si, ■ ■ ■ ,S n given real constants, 
with oi, • • • , a n given positive constants, with suitable arrays (qjf)i<i,j<n of real constants, the 
"descending order statistics" notation 

max Xj(t) =: X {1) (t) > X (2) (t) > ■ ■ ■ > X (n _ 1} (t) > X (n) (t) := min Xj(t) , 

and lexicographic breaking of ties. This system (1.8) exhibits both features of rank-dependent 
characteristics and skew-elastic collisions that are manifest in (1.2), (1.3), but involves several 
particles rather than just two. 

2 On Semimartingale Local Time 

Let us recall the notion of a continuous, real-valued semimartingale 

X(-) = X(0) + M(-) + C(-), (2.1) 
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where M(-) is a continuous local martingale and C(-) a continuous process of finite first variation 
such that M(0) = C(0) = . The /ora/ rime f ) accumulated at a given "site" f e R over 

the time-interval [0, i] by this process, is 

1 /"* 

L x (t;0 := lim — / l { ^ x(s)<+£} d(X)(s) 

£-10 Z £ Jg 

= (X(t)-£) + -(X(0)-£) + - f l Ws) >£ } dX( S ), (2.2) 

Jo 

where = (M)(-) . For every fixed £ G R this defines a nondecreasing, continuous and 

adapted process L x (- ; £) which is flat off the set {t > : X(t) = £} , namely 

POO POO 

/ l{x(t)^} dL x (t; £) = 0; and we have also the property / l{x(t)=£} d(X)(£) = 0. 
Jo Jo 

(2.3) 

On the other hand, for each fixed T e (0, oo) the mapping £ i->- L X (T; £) is almost surely 
RCLL (Right-Continuous on [0, oo), with Limits from the Left on (0, oo)) paths, and has jumps of 



size 

fT 



L x (T;0-L X (T;£-)= [ l {x(t)= ^ } dX(t) = [ l {x(t)=e} dC(f) . (2.4) 
We shall employ also the notation 

L X (T; := j (L X (T; £) + L X (T; £-)) (2.5) 

for the so-called "symmetric local time" accumulated at the site £ over the time interval [0, T] . 
For these local times we prefer to use the simpler notation 

L x {-) = L x (.;0), L x {.) = L x (.;0-), L x (-) 0) (2.6) 

when we evaluate them at the origin £ = , and note 

L x (-) = L- x (-) = lim -L f l {0 > x{t)> _ £} d(X)(t). (2.7) 

eiO Z £ J q 

Finally, we recall the occupation time density formulae 

/ h(X(t))d(X)(t) = 2 / L x (-;0H0^ = 2 / L x (-;0H0^, (2.8) 
Jo Jr Jr 

valid for every Borel measurable function h : R — >■ [0, oo) , as well as the ITO-TANAKA formulae 

/(*(•)) = /(X(0))+ / ZT/(X(t))dX(t) + / ( 2 -9) 

Jo Jr 

/(*(•)) = /(*(0)) + i / (^ + /(X(t)) + D-/(X(t)))dX(t)+ / L*(-;0/W (2-10) 

z Jo Jr 

Here / : R — )■ R is the difference of two convex functions, D ± f(-) denote its derivatives from 

left and right, and /"(•) denotes its second derivative measure. 
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2.1 Tanaka Formulae 

For a continuous, real-valued semimartingale X(-) as in (2.1), and with the conventions 

sp(x) := l( ,oo)(aj) - l(-oo,o)(aO , sgn (x) := l(o,oo)(z) - l(_oo,o](aO , x G R (2.11) 

for the symmetric and the left-continuous versions of the signum function, we obtain from (2.9), 
(2.10) the Tanaka formulae 

\X{-)-i\ = \X(0)-£\+ [ sgn(Jf(t)-£)dJf(t) + 2L x (-;0 (2.12) 

Jo 

= |X(0)-e|+ / W(X(t)-£)dX(t) + 2L x (-;0- (2-13) 
Jo 

Applying (2.12) with £ = to the continuous, nonnegative semimartingale |X(-)| , then compar- 
ing with the expression of (2.12) itself, we obtain the companion 

2L x (-)-L\ x \-) = [ l {x(t)=0} dX(t) = [ l {x{t)=0} dC(t) (2.14) 
Jo Jo 

of the property (2.4). For the theory that undergirds these results we refer, for instance, to KARATZAS 
& Shreve Karatzas, I. & Shreve, S.E. (1991) (1991), section 3.7. 

3 Analysis 

Let us suppose that such a probability space as stipulated in section 1 has been constructed, and 
on it a pair -Bi(-) , B 2 (-) of independent standard Brownian motions, as well as two continuous, 
nonnegative semimartingales -Xi(-), X 2 (-) such that the dynamics (1.2)-(1.3) are satisfied. We 
import the notation of Fernholz et al. Fernholz, E.R., Ichiba, T., Karatzas, I. & 
PROKAJ, V. (2012) (2012) : in addition to (1.1), we set 

v = g — h , y = xi — X2 , z := X\ + x% > , r\ = X\ V x% , r 2 = X\ A x 2 , (3.1) 

and introduce the difference and the sum of the two component processes, namely 

Y(-) := Xi(-) - X 2 (.) , Z(-) := X 1 (-) + X 2 (-) . (3.2) 

3.1 Auxiliary Brownian Motions 

We introduce also the two planar Brownian motions (Wi(-), W 2 {-)) and (Vi(-), V 2 (-)), given 
respectively by 

W x {-) := / l {Y(t)>0} dB 1 (t)- [ l {Ymo} dB 2 (t), (3.3) 
Jo Jo 

W 2 (-) := [ l {Y{t) < 0} dB 1 (t)- [ l {Y{t)>0} dB 2 (t) (3.4) 
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and 



Vi(-) := / l {Y (t)>o}dB 1 (t)+ [ l {Ymo} dB 2 (t), (3.5) 
Jo Jo 

V 2 (-) := / l {Ymo} dB 1 (t)+ [ l {Y{t)>0} dB 2 (t). (3.6) 
Jo Jo 

Finally, we construct the Brownian motions W(-), V(-), Q(-) and W^^), V b (-), U\-) as 
W(-) := pW 1 (-) +aW 2 (-), V(-) := P V X {-) + aV 2 {-) , Q(-) := a V x {-) + p V 2 {-) , (3.7) 

W\-) :^ P W 1 {-)-aW 2 {-), V\-) := P V 1 (-)-aV 2 (.) , U\-) := a W 1 (-)- p W 2 (-) ; (3.8) 

we note the independence of Q(-) and W(-) , the independence of Q(-) and V (•) , and observe 
the intertwinements 

VS-(-) = [ sgn(Y(t))dWj(t) (j = 1,2), = / sgn(y(t))dW(t) (3.9) 

Jo Jo 



and 



V(-)=/ sgn(F(t))diy b (t), Q(-)=/ sgn(r(t))df/ b (t). (3.10) 
Jo Jo 



3.2 The Difference and the Sum 

After this preparation, we observe that the difference Y(-) and the sum Z(-) from (3.2) satisfy, 
respectively, the stochastic integral equation 

Y{-) = y-X [ sgn(F(t))dt+(l-C)^ y (-)-(l-^)^ y (-) + ^(-) (3.H) 
Jo 

which involves both the right- and the left- local time at the origin of its solution process Y(-) , 
and the identity 

Z(t) = z + ut + V(t) + (l-()L Y (t) + (1 -r/) L y (t) , 0<t<oo. (3.12) 
We have used here the notation of (1.1), (3.1), (1.4), as well as 

We note from (2.3) and (3.11) that Y(-) is a continuous semimartingale with 

/*oo POO 

/ l {Y(t)=0} dt = l {m=0} d{Y){t) = 0, a.s., (3.14) 
Jo Jo 

and that on the strength of (2.4), (2.3) we have 

L Y (- ; 0) - L Y (- ; 0-) = / l {y(t)=0} ((1 - C) dL Y (t; 0) - (1 - V ) dL Y (t; 0-)) 

Jo 



= (l-C)L y (-;0)-(l-r/)L y (-;0-) 

or equivalently 

CL Y {.) = V L Y (.). (3.15) 
From this relationship and (2.3)-(2.5), (2.14) we obtain 

2L Y (-) = and L Y '(•) = aL |y| (-), L Y {-) = (1 - a) L |y| (-) , (3.16) 

where we introduce as in (1.5) the "skewness parameter" 

a:=^—. (3.17) 

v + C 

• With this notation, and recalling (3.14) and (3.16), we see that equation (3.11) takes the form 

Y(-) = y-\[ sp(y(t))dt + 2(2a-l)L y (-) +W{-) (3.18) 
Jo 

of the equation for a Skew Brownian Motion with Bang-Bang drift (Skew Bang-Bang Brownian 
Motion, or SBBBM for short). This is a very close relative of the Skew Brownian motion, that 
was introduced by Ito & McKean Ito, K. & Mc Kean, H.P., Jr. (1963) (1963), Ito, K. 
& McKean, H.P., Jr. (1974) (1974) and was furhter studied by Walsh Walsh, J.B. (1978) 
(1978), Harrison & Shepp Harrison, J.M. & Shepp, L.A. (1981) (1981) ; see Lejay 
Lejay, A. (2006) (2006) for a comprehensive survey. 

The diffusion process Y(-) of (3.18) is studied in detail in section 6. It is a strong Markov 
and Feller process, whose transition probabilities can be computed explicitly; see (6.13)-(6.15) 
below. In particular, it is shown in section 6 that, for < a < 1 , the stochastic equation (3.18) 
has a pathwise unique, strong solution and that the filtration identities 

d Y (t)=$ w (t), ViG[0,oo) (3.19) 

hold. Here and in what follows, given a process H : [0, oo) x fi — > R. d with values in some 
Euclidean space and RCLL paths, we shall use the convention 

F H = {S H (t)}o<t<oo 

for the smallest filtration to which H(-) is adapted that satisfies the "usual conditions" of right 
continuity and augmentation by sets of measure zero. 

• Similarly, and with the notation of (3.17), the expression (3.12) takes the form 

X x (t) +X 2 (t) = Z(t) = z + ut + V(t) + 2 (1 -0) L Y (t), 0<t<oo (3.20) 
where, as in subsection 1.1, we set 

£:=^i^. (3.21) 
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Remark 3.1. It is clear from (3.15) that the equation (3.11) can be written as 

Y(-) = y - X f sgn(F(t)) dt + L Y {-) - L Y {-) + W(-) . (3.22) 
Jo 

To wit: skew Brownian motion with bang-bang drift solves the equation (3.22), for any value 
a E [0, 1] of its skewness parameter. We conclude that uniqueness in distribution fails for this 
equation (3.22), thus also for the equation (3.11) that governs the difference Y(-) = Xi(-) — X 2 (-) 
when r] = ( = . 

In particular, uniqueness in distribution cannot possibly hold for the system (1.2), (1.3) in this 
case rj — C = . 

Remark 3.2. When rj = -( ^ we get L Y '(•) + L Y (-) = from (3.15), thus 

L Y (.) = L Y (.) = 0, (3.23) 
and the equation (3.11) takes the form of Brownian motion with bang-bang drift 

Y(-) = y - X f sgn(F(t)) dt + W(-) . (3.24) 
Jo 

This diffusion process was studied in detail by Karatzas & Shreve Karatzas, I. & Shreve, 
S.E. (1984) (1984), who computed its transition probabilities and the joint distribution of the triple 
( Y(t), L Y (t), J ' l{y( s ) >0 } ds) . This diffusion does accumulate local time at the origin: indeed, 
on the strength of (2.4), (2.3), we have 

L Y (.)-L Y (-)= [ l {Y{t)=0} dY(t) = X [ l {Y(t)=0} dt = X [ l {Y(t)=0} d(W)(t) = 
Jo Jo Jo 

almost surely, but also F(L Y (t) = L Y (t) > 0) > for every t E (0, oo) ; this contradicts (3.23). 
In fact, we have F(L Y (t) = L Y (t) > 0) = 1 for y = . 

We conclude that the equation (3.11) for the difference Y(-) = Xi(-) — Xa(-) has no solution 
in the case i] = — £ ^ . Thus, the system (1.2), (1.3) cannot possibly have a solution in this case. 

3.3 Auxiliary Systems 

From the equations of (3.11), (3.12) and using the notation in (3.2)-(3.7), we obtain a system of 
stochastic differential equations 

dJfi(t) = ^l{Xi(t)<x 2 (t)} - hl {Xl {t)>x 2 (t)} S j dt + pl {Xl (t)>x 2 {t)}dWi(t) 

+ a l {Xl (t)<x 2m dW 2 (t) + dL M2 (t) + dL x ^(t) , (3.25) 

dX 2 (t) = (g l{Xi(t)>x 2 (*)} - h l{xi(t)<x 2 (t)}) dt - pl {Xl ( t )<x 2 (t)} dW^t) 

1 C 1 
- a l {Xl{ t)>x 2{ t) } dW 2 (t) + dL x ^ x '\t) + dL x ^(t) , (3.26) 
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quite similar to that of (1.2), (1.3), but now driven by the planar Brownian motion (Wi(-), W 2 (-)) . 
In a totally analogous manner, we obtain also the system 

dXi(t) = r^l{xi(t)<x 2 (t)} - hl{x 1 {t)>x 2 {t)}j dt + pl{x 1 (t)>x 2 (t)}dV 1 (t) 

1 C 1 
+ a l {Xl{ t)<x 2 (t)} dV 2 (t) + dL x ^{t) + dL x ^(t) , (3.27) 

dX 2 (t) = (g l{xi(t)>x 2 (t)} - frl{Jfi(t)<x 2 (t)}) dt + pl{ Xl (t)<X2(t)}dVi(t) 

+ a l { x l{ t)>x 2m dV 2 (t) + dL*-*(t) + dL*~*(*) , (3-28) 

now driven by the planar Brownian motion (Vi(-), V 2 (-)) . 

3.4 Skew Representations 

In light of the Tanaka formula (2.13), of the equation (3.18) for the semimartingale F(-) , and of 
the last intertwinement in (3.9), we represent now the size of the "gap" between Xi(t) and X 2 (t) 
as 

\Y(t)\ = \y\ - \t + V\t) + 2L Y (t) (3.29) 
= \y\-Xt + V\t) + L lYl {t), 0<t<oo. 

With the help of (3.9), (3.16), let us write the first Brownian motion in (3.8) as 

W\-) = 7 W(-) + 5U\-) , where 7 := p 2 - a 2 , 5 := ^1 - 7 2 = 2pa. 
With this notation, and with the help of (3.10), the Brownian motion V(-) in (3.7) takes the form 

V(t) = jV\t) + SQ(t) = -f(\Y(t)\-\y\ + \t-2L Y (t))+5Q(t), < t < 00 . 

We recall here from (3.10) the standard Brownian motion Q ( • ) which, being independent of W( • ) , 
is also independent of the process Y(-) in light of (3.19). 

In conjunction with Xi(t) — X 2 (t) = Y(t) and the representation (3.20) for Xi(t) + X 2 (t), 
and with the notation 

H := — (z/ + A 7 ) = gp 2 -ha 2 , 
we obtain from this expression the skew representations for the component processes themselves 
X 1 (t)=x 1 +pt + p 2 {Y + (t) - y + ) - a 2 (Y-(t) - y~) + (l - /? - 7) L Y (t) + paQ(t) (3.30) 
X 2 (t) = x 2 + pt - a 2 (Y + (t) - y + ) + p 2 - y~) + (l - (3 - 7) L y (t) + p<rQ(t) (3.31) 

in terms of the paths of the skew Brownian motion process K(-) with bang -bang drift, and of the 
independent Brownian motion Q(-) . In particular, this shows that uniqueness in distribution holds 
for the system of stochastic differential equations (1.2), (1.3). 

Similar reasoning shows that uniqueness in distribution holds also for each of the systems 
(3.25), (3.26) and (3.27), (3.28). 
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Remark 3.3. It is clear from (3.29) that the absolute value of the skew Brownian motion with 
bang-bang drift in (3.18), for any value a 6 [0, 1] of the skewness parameter, is Brownian motion 
with drift -A and reflection at the origin. Arguing as in Walsh Walsh, J.B. (1978) (1978) , 
Proposition 1, one can conclude that every diffusion process Y(-) , for which \Y(-)\ is Brownian 
motion with drift —A and reflected at the origin, is a skew Brownian motion with bang-bang drift. 

4 Synthesis 

We reverse now the steps of the analysis in section 3. Let us start with a filtered probability space 
P) , F = {$(t)}o<t <0 o and with two independent, standard Brownian motion Wi(-) , W 2 (-) 
on it; we shall assume F = F ( Wl > w ?) , i.e., that the F is the smallest filtration satisfying the usual 
conditions, to which the planar Brownian motion (Wi(-), W 2 (-)) is adapted. 

With given real constants Ci ? C2 , Vi > V2 and nonnegative constants g , h, p, a that satisfy 
(1.1) and (1.5), with a given vector (x±, x 2 ) G , and with the notation of (3.1), we construct the 
pairs of independent Brownian motions 

W(-) : = pWx(-) + aWz(-), U\-) : = aW x {-) - P W 2 {-) (4.1) 

and 

[/(-) : = aW x {-) + P W 2 (-), W\-) : = pWi(-) - aW 2 (-) (4.2) 

as in (3.8), (3.7). Clearly, F^uW,) = F (w,u») = F (u,w») 

We construct also the pathwise unique, strong solution Y(-) of the stochastic equation (3.18) 
driven by the Brownian motion W(-) introduced in (4.1). With the process Y(-) thus in place, 
we introduce by analogy with (3.9) the independent Brownian motions 

V 1 (-)= [ sgu(y(t))dWi(t), V 2 (t) = - [ sga(Y(t))dW 2 (t), (4.3) 
Jo Jo 

and by analogy with (3.7), (3.8) the two additional pairs of independent Brownian motions 

V(-) := pVxO) + aV 2 (-), Q\-) := aVx(-) - P V 2 (-) (4.4) 

and 

Q(-) := aVx{-) + pV 2 (-), V\-) : = pV x {-) - aV 2 {-) . (4.5) 

• We introduce also the continuous martingales 

Mx(-) := J (pl { r W > }dWi(t) + (7l { y (t) <o } dW 2 (t)) (4.6) 

(pl {Y (t) >0 } dVx(t) + crl{ Y (t)<o} dV 2 (t)j , 




12 



M 2 (-) := - J (p l { y W <o} dWi(t) + a l { y (t) >o} dW 2 (t)) (4.7) 

= y (p l{Y(t)<o} dVi(t) + cr l{y(t)>o} dV£(t)) , 
with (Mi,M 2 ) (-) = and quadratic variations 

(Mi)(0= / (p 2 l { y W >o} + ^ 2 l{y( t )<o})dt, (M 2 )(0=/ (p 2 l { y (t )<o} + ff 2 l { y( t )>o})dt. 
Jo Jo 

There exist then independent Brownian motions £>i(0 , -B 2 (0 on our filtered probability space 
(Q, P) , F = {S r (t)}o<t<oo j so the continuous martingales of (4.6), (4.7) are cast in their DOOB 
representations as 

M x {-) = J (pl{y( ( )>o}+ffl{r(()<o}) dB x {t) , M 2 (0 = J (pl{y(()<o}+^l{r(t)>o}) dB 2 (t) 

(4.8) 

in terms of independent Brownian motions B x (-) , B 2 (-) ; for instance, by taking 

#i(0 = J ( l{r(i)>o} dW x (t) + l { y (t) <o} dW 2 (t)) , (4.9) 
B 2 {-) = -J (}{Y(s)<o} dW x {t) + 1{y W >o} dW 2 (t)) . (4.10) 

• Finally, we introduce the continuous, F— adapted processes 

X 1 (-):=x 1 + j (gl {Ymo} -hl {Y{t)>0} )dt + M x {.) + L Y {■) + ^—^L Y _{-) 

(4.11) 

and 



1-^2 y ,. 



X 2 (-):=x 2 + J (gl {Y(t)>0} -hl {Ymo} )dt + M 2 (0 +^^L Y (-)+ ^ 

(4.12) 

It is now easy to check Xi(0~ X 2 (0 = Y(-) ; and from this, that the vector process (Xi(0, X 2 (-)) 
solves the system (1.2)-(1.3), as well as the systems (3.25)-(3.26), (3.27)-(3.28). It is also straight- 
forward to verify the skew representations of (3.30), (3.31). 

Remark 4.1. We note that the vector process (Xl(0, -^2(0) solves also the system of stochastic 
equations 

dX x (t) = l {Xim x 2m (gdt + crdB x (t)) (4.13) 
+ l {Xl (t)>x 2 (t)}{- hdt + pdB^t)) +K 1 dL^- x ^(t), 

dX 2 (t) = l {Xl{t)>X2{t)y (gdt + adB 2 {t)) (4.14) 
+ l{ Xl (t)<x 2 (t)} {-hdt + p dB 2 (t)) + k 2 dL^~ x ^(t) , 
where 2 Kj := a (1 — Q) + (1 — a) (1 — rjj) , j = 1, 2 or equivalently 

ki = a-(p/2), k 2 = 1 -a- (£/2) . (4.15) 
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4.1 Ranks 



Let us introduce explicitly the ranked versions (leader and laggard, respectively) 

i2 x (.) := Xi(-) V X 2 (-) , := X x (-) A X 2 (-) (4.16) 

of the components of the vector process (Xi(-), X 2 (-)) constructed in (4.11), (4.12). From (3.12) 
and (3.29), it is rather clear that we have 

Ri(t) + R 2 {t) = Xi(t)+X 2 (t) = ri + r 2 + ut + V(t) + (1 -/3)L |y| (t), < t < oo 

R 1 (t)-R 2 (t) = \X!(t) - X 2 (t)\ = \Y(t)\ = \y\-\t + V\t) + L^(t), (4.17) 

and these representations lead to the expressions 

R^t) = n - ht + pVi(t) + (1 - (13/2)) L Rl ~ R2 (t) , 0<t<oo (4.18) 

R2(t) = r 2 + gt + aV 2 (t)-(p/2)L R '- R *(t), < t < oo . (4.19) 

A few remarks are in order. The equations (4.18), (4.19) identify the processes Vi(-) and 
V 2 (-) of (3.5), (3.6) as the independent Brownian motions associated with the diffusive motion 
of the ranked particles, the "leader" R\(-) and the "laggard" R 2 (-) , respectively; whereas the 
independent Brownian motions in (1.2) and B 2 (-) in (1.3) are associated with the specific 

"names" (indices, or identities) of the individual particles. On the other hand, with the help of 
the theory of the Skorokhod reflection problem (e.g., Karatzas & Shreve Karatzas, I. 
& Shreve, S.E. (1991) (1991) , page 210), we obtain from (4.17), (2.3) the identification of the 
"collision local time" 

LRi-R*( t ) = L |y| (t) = max (- \y\ +Xs-V\s)) + , < t < oo . (4.20) 

0<s<t V / 

Let us also observe that, in the non-degenerate case pa > , the equations (4.17)-(4.20) and the 
second equation in (4.5) give the filtration comparisons 

#(W(t) = & R M(t) , < t < oo , (4.21) 

3 V \t) = d lYl (t) c 3 rV( t ) | <t<oo, (4.22) 

where the inclusion is strict, due to the fact that the process Y(-) changes its sign with positive 
probability during any time-interval [0, t] with t > . 

4.2 Filtration Comparisons, Weak and Strong Solutions 

We have the following straightforward analogues of Propositions 4.1, 4.2 and of Theorems 4.1, 
4.2 in Fernholz et al. Fernholz, E.R., Ichiba, T., Karatzas, I. & Prokaj, V. (2012) 
(2012). 
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Proposition 4.1. In the degenerate case a = 0, thus p = 1 in light of (1.1 ), we have the relations 

£(*i,-Ra)( t ) = $V(q = g\x,-x 2 \^ c ^^i-^(t) = j^(t) = £(*i.*0(t) (4.23) 

/or every < t < oo , where the inclusion is strict. 

In the special case (3 = 1 of (5.2) we have in addition cr(V(t)) = cr(Xi(t) + X 2 (t)), thus 
also $ v (t) = $ Xl+X2 (t) , for every < t < oo. 

Proposition 4.2. In the non-degenerate case pa > 0, we have for every < t < oo the filtration 
relations 

#M)(t) = $(*A)(t) = = 5 (|y| ' Q) (t) (4.24) 

cgr(^)(t) = j( y ^(t) = $WM(t) = d (Xl ' X2 \t), 

where the inclusion is strict. 

Theorem 4.1. The system of stochastic differential equations (1.2), (1.3) is well-posed, that is, has 
a weak solution which is unique in the sense of the probability distribution. The same is true for 
each of the systems of equations (3.25), (3.26) and (3.27), (3.28). 

On the other hand, the system of stochastic differential equations (3.25), (3.26) admits a strong 
solution, which is therefore pathwise unique; whereas the system (3.27), (3.28) admits no strong 
solution. 

Theorem 4.2. The system of stochastic differential equations (1.2), (1.3) admits a pathwise unique, 
strong solution; in particular, the filtration identity $( Bl > B2 >(t) = $( Xl > X2 \t) holds for all < 
t < oo . 

Likewise, the system of equations (4.13), (4.14) has a pathwise unique, strong solution. 

Proof. Repeating almost verbatim the arguments in the proof of Theorem 5.1 in Fernholz ET 
al. Fernholz, E.R., Ichiba, T., Karatzas, I. & Prokaj, V. (2012) (2012), the question 
boils down to whether the filtration comparison 

$ Y (t) C S^ 1 '^) , V < t < oo (4.25) 

holds. To decide this issue, we write the equation (3.18) as driven by the pair B 2 (-)) ; in 

other words, we use (3.7) and (3.3), (3.4) to express the skew Brownian motion F(-) with bang- 
bang drift as solution of a stochastic differential equation driven by the planar Brownian motion 
(Bi(-), B 2 (-)) . Since this equation does admit a weak solution which is unique in distribution, the 
issue is whether this solution is also strong, that is, whether (4.25) holds. 

This question is easy to settle in the isotropic case p = a = l/y/2; then W(-) = (-Bi(-) — 
B 2 {-)) I \/2 , and the comparison (4.25) follows from the strong solvability of the equation (3.18) 
proved in section 6: $ Y {t) = $ w (t) C ^^' B2 \t) holds for all < t < oo , by virtue of (3.19). 

In the non-isotropic case p ^ o , we write (3.18) as the extended skew Tanaka equation 

Y(t) = y+ [ t m(Y(s))dp(s)-P±^V(t)+2(2a-l)L Y (t), 0<t<T 

V 2 Jo V 2 

(4.26) 
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with T G (0,oo) arbitrary but fixed. Here 

/3i(0+/3 2 ( 



/3(- 



0(0 



are given standard, independent Brownian motions after an equivalent change of probability mea- 
sure, and 

P&) := Btf)--^-, 0<t<T (2 = 1,2). 
p-o 

• Let us suppose that Fi(-) and Y 2 (-) are two solutions of the equation (4.26), defined on the 
same probability space and with respect to the same, independent standard Brownian motions 
Bi(-), B 2 {-) . Following Le Gall (1983), we shall show L Yl ~ Y2 {-) = ; we shall then argue that 
this implies also Y\(-) = Y 2 (-) , a.s. 

To this end, we consider the difference D(-) := Yi(-) — Y 2 (-) , as well as the linear combina- 
tions 

Z (u) (-) := (1 -u)Y x {-) + uY 2 (-) for 0<it<l; 

we introduce also a sequence {fk}ken C C X (]R) of continuous and continuously differentiable 
functions that converge to /oo(-) : = sgn(-) pointwise, and satisfy sup fc6N ||/fe|| ry < 00. Since 
limsup fc ||/fe|| T y < || /oollyy obviously holds, this is only possible if /oo(-) is of bounded variation, 
and in this case an approximating sequence is easily obtained, e.g., by mollifiers. As in the proof 
of Theorem 8.1 of Fernholz et al. Fernholz, E.R., Ichiba, T, Karatzas, I. & Prokaj, 
V. (2012) (2012), for every 5 > , T > , k > 1 ,we establish then 



E 







IfkjY^s)) - f k (Y 2 (s))\ 
Y 1 ( 8 )-Y 2 (s) 



l{yi( s )-y 2 (s)>5} dt 



< c a ||iy Ty -supE(2L»(T,0) 



here L^T, £) is the symmetric local time of Z^ u \-) accumulated at the site £ G R over the 
time interval [0, T] , and c x is a constant chosen independently of k,u,5 . Letting k t 00 and 
5 I , we obtain 



E 



D(t) 



l {D{t>0} d(D)(t) 



< 2E 









Y(t)-Y 2 (t) 1{m>0} dt 



< 2 Cl \\f oo \\ TV -snpE(2L^(T,0 



Now the Cauchy-Schwartz inequality, the ITO isometry, and the TANAKA formula (2.13) ap- 
plied to Z^ u \-) , allow us to estimate 

E(2L (u) (T;0) < E\Z {u) (T) - Z (u) (0)| + \E((Z {u) )(T))] 1/2 

+ 2(2a-l)(uE(L Yl (T)) + (1 - u)E(L Y *(T))) 



< 2 



[E((Z (u) )(T))] 1/2 + 2 (2a - 1) (uE(L Yi (T)) 



u 



E(L Y \T)) 
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The last term is bounded uniformly in (£,«), since (Z^)(t) < c 2 t and E(L y '(T)) < C3 , for 
% — 1, 2 and for some constants C2 , C3 that do not depend on (£, u) . Thus, we obtain 



E 



T 



1 



D(t) 



l {D(t)>0} d(D)(t) 



< 00 , < T < 00 



(4.27) 



Using Lemma 1.0 of Le Gall LeGall, J.F. (1983) (1983) (see also Exercise 3.7.12, pages 
225-226 in Karatzas & Shreve Karatzas, I. & Shreve, S.E. (1991) (1991)), we verify that 
(4.27) gives L D (-) = . By exchanging the roles of li(-) and Y 2 (-) , we obtain also L~ D (-) = 
L Y2 ~ Yl (-) = , as well as L D {-) = . Furthermore, we note that on the strength of Corollary 2.6 
of OUKNINE & Rutkowski OUKNINE, Y. & RUTKOWSKI, M. (1995) (1995) this implies that 
the symmetric local time L M {-) of the maximum 

M(-) := ¥,{■) V Y 2 (.) = Yl(.) + (Y 2 Q - F 1 (-)) + 

is given as 

L M (.) :=L Y ^(-)= [ l {Y2{t)<0} dL Y (t)+ [ l {Yl{t)<0} dL Y *(t). 

Jo Jo 

We combine now these results with the TANAKA formula, to obtain the dynamics of the maximum 

M(-)=y+[ l {Yl{t)>Y2m dY 1 (t)+ [ l {Yl{t)<Y2m dY 2 (t) + L Y ^(-) 
Jo Jo 

= y + I sp(M(t)) d(3(t) - £±1 + 2 (2q - 1) L M {.) , 

and observe that these are the same as those of (4.26). But uniqueness in distribution holds for the 
equation (4.26), so the distribution of the process M(-) is the same as that of Yi(-) ; and of course 
we have M(-) > Y x {-) a.s. This implies M(-) = Yi(-) , thus Yi(-) = Y 2 (-) a.s. 

Therefore, the solution to (4.26) is pathwise unique, hence also strong by the theory of Ya- 
mada & Watanabe (e.g., Karatzas & Shreve Karatzas, I. & Shreve, S.E. (1991) 
(1991), pages 308-311). □ 



5 Some Special Cases 

When a = 1/2 , that is, r] = ( ^ or equivalently 

771-772 = C2-C1 7^ 2, (5.1) 

the equation (3.18) for the difference Y(-) = Xi(-) — X 2 (-) becomes that of Brownian motion 
with bang-bang drift 

Y(t)=y-\f sgn(Y(s)) ds + W{t) , < t < 00 
Jo 
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as in (3.24). In this special case jj = ( ^ and with a = p, the existence and uniqueness of 
(1.2)-(1.3) can be shown also by direct application of Theorem 3.5 of Sznitman & Varadhan 
Sznitman, A.S. & Varadhan, S.R.S. (1986) and a Girsanov's change-of-measure, with 
the aid of the local time relationships (3.16). 

On the other hand, when f3 — 1 or equivalently 

V{l-C) = C(l-»7), (5-2) 

we observe from (3.20) that the sum Xi(-) + X 2 (-) is just standard Brownian motion with drift 
v = g-h. 

Let us single out now, and study, some more interesting special cases. 

5.1 Perfect Reflection for Individual Particles Upon Collision 

Suppose a = 1 , or equivalently £ = and 77 ^ from (1.5), that is 

C2-C1 = 2 ^ 771-%, (5.3) 

and that x 1 > x 2 . We see then from (3.16) that we have L Y (-) = L Xi ~ Xl (-) = , and that the 
equation (3.18) becomes 

Y(t) = y-Xt + W(t) + L Y (t) } 0<t<oo. 

From the theory of the Skorokhod problem (e.g., Karatzas & Shreve Karatzas, I. & 
Shreve, S.E. (1991) (1991), pages 209-210) we conclude 

L Y (t) = max (-y + Xs-W(s)) + , < t < 00 

0<s<t v ' 

thus Y(-) > 0, and that strength and pathwise uniqueness hold; for more general results along 
these lines see Chitashvili & Lazrieva Chitashvili, R.J. & Lazrieva, N.L. (1981) (1981). 
It is also clear from the last two displayed equations, that the filtration identity # y (t) = $ w (t) , < 
t < 00 in (3.19) also holds. 

In this case, then, when the particles collide, the trajectory Xi(-) of the first particle bounces 
off the trajectory X 2 (-) of the second particle as if this latter were a perfectly reflecting lower 
boundary. We can visualize the situation by saying that, under the conditions of (1.5) and (5.3), 
the second particle is "heavy" (unaffected by collisions), whereas the first particle is "light" in that 
it bounces off (reflects perfectly) when colliding with the heavy particle. 

• The "symmetric" situation obtains for a = , that is ( ^ and rj = or equivalently 

C2-C1 + 2 = Tfc-Tfc ; (5.4) 

in this case and again with x\ > x 2 , when the two particles collide, the second particle bounces 
off the first as if this latter were a perfectly reflecting upper boundary; it is the first particle that is 
now "heavy", and the second that is "light". 
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5.2 Frictionless Collision 

It follows also from (3.16) that the local times disappear entirely in (1.2) when we have the con- 
figuration of parameters (1 — £i)a + (1 — 771) (1 — a) = , or equivalently 

(l-Ci)r/ + (l-r/ 1 )C = 0; (5.5) 

in this case the trajectory of the first particle crosses that of the second without "feeling it", that is, 
without being subjected to any local time drag. 

Similarly, the second particle crosses the first in the same frictionless manner, that is, the local 
times disappear entirely in (1.3), if 

(1-C 2 )?7 +(1 - 77 2 K = 0. (5.6) 

• If both (5.5) and (5.6) hold, then all such crossings are completely frictionless. We note that 
(5.5) and (5.6) are both satisfied, if and only if 

m + Ci = m + C2 = 2 (5.7) 

holds. This condition implies 77 = ( (so when this common value is nonzero we are in the 
case a = 1/2 mentioned at the start of the section), and is obviously satisfied in the special case 
771 = d = r/ 2 = C2 = 1 studied by Fernhoz et al. Fernholz, E.R., Ichiba, T., Karatzas, 
I. & Prokaj, V. (2012) (2012). However, (5.7) holds also for other configurations of parameters, 
for instance £1 = 772 = 1/2 , 771 = (2 = 3/2 . 

The condition (5.7) gives the value (3 = 1 for the parameter of (3.21); back in (4.18), (4.19), 
this implies that the collision local time L Rl ~ R2 (-) "gets apportioned equally to the ranks". 

5.3 Elastic Collisions 

Beyond these two extremes of perfect reflection and frictionless collision - that is, for all other 
configurations of parameters - we have collisions that are "elastic": neither completely frictionless, 
nor perfectly reflecting. 

5.4 Brownian motion reflected on an independent Brownian motion 

Finally, let us consider the case (3 = or equivalently 77 £ + £ 77 = , that is 

2 (Ci + C2 + V1 + V2) = (Ci + Ca) (m - V2) - (7/1 + 772) (Ci - C2) (5.8) 

in light of (3.21) and (1.4), (3.13). This happens, for instance, when (1 = 3/4, ( 2 = 9/4, 
771 = —4/3 , 772 = —8/3 ; in this case we have a = 4/7 and of course j3 = . 

Under the condition (5.8), the laggard in (4.19) feels no pressure (local time drag) from the 
leader; it just evolves like Brownian motion with variance a 2 and nonnegative drift. On the other 
hand, the leader in (4.18) evolves like an independent Brownian motion with variance p 2 and 
nonpositive drift, reflected off the trajectory of the laggard. Such a process has been studied by 
Burdzy & Nualart Burdzy, K. & Nualart, D. (2002) (2002) (see also SOUCALIUC ET 

AL. SOUCALIUC, F., TOTH, B. & WERNER, W. (2000) (2000), SOUCALIUC & WERNER 
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Soucaliuc, F. & Werner, W. (2002) (2002)); here it arises as a special case of the ranked 
system (4.18), (4.19) for the particles whose motions are governed by the equations (1.2), (1.3). 

We have in this case (3 = an interesting fusion: the "perfect reflection" we saw in subsection 
5.1, and the "frictionless motion" of subsection 5.2, are occurring here simultaneously - not for the 
motions of the individual particles, however, but rather for the motions of their ranked versions, 
the leader Ri(-) and the laggard -R2O) , respectively. To put it a little differently: starting with two 
particles that undergo skew-elastic collisions one is able, under the conditions of (1.5) and (5.8), 
to "simulate a heavy particle" (the laggard) and a "light" particle (the leader). 

• The "reverse" situation obtains when (3 = 2 or equivalently + (,rj = 2 (?? + C) > that is 

2(Ci + C 2 + ^i + r/ 2 ) = 4 (4 + Ci-C 2 -?7i+?72) + (C1 + C2) (771-^2) -(t7i + %)(Ci-C 2 ) ! (5.9) 

then it is the trajectory of the laggard (now the "light" particle) that gets reflected off that of the 
leader (now the "heavy" particle). This happens, for instance, when £1 = 3/2 , £2 = 3 , r)i = 7/3 , 
772 = 1 ; in this case we have a = 4/7 and (3 = 2. 

5.5 Some Simulations 

The pictures (Figures 1-4) that follow present simulations of the processes Xi(t) (in black) and 
X 2 (t) (in red) for t 6 [0, 1] , in black and red, respectively, with drifts g = h = 1 in the 
degenerate case p = . 




Figure 1: Ci = C2 = Vi = V2 = 1 ; a = 1/2 , /? = 1 . 
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Figure 4: Ci = 0, C2 = 2, 771 = 772 = !; « = 1, /3 = 1. 



6 Skew Brownian Motion with Bang-Bang Drift 

We study here the stochastic differential equation (3.18) for the skew Brownian motion with bang- 
bang drift 

b(y) = -Asgn(y), y eR (6.1) 

for some given constant A > 0, with the notation of (2.11), and with skewness parameter a G 
[0, 1] . The cases a = and a = 1 have been discussed already in subsection 5.1, so we focus 
here on the range < a < 1 . 

For this range of values of the skewness parameter, we choose to write the equation (3.18) in 
terms of the right-continuous local time of the unknown process at the origin, namely 

Y(-) = y - A / sgn(Y(t)) d* + W(-) + ^— ^ L Y {-) . (6.2) 
Jo " 

This equation is of the more general form 

Y(-) = 7/o+ f r{Y{t))dW(t)+ [ L Y (-,0"(dO, (6.3) 
Jo Jr 

with dispersion r(y) = 1 and measure 

2 a — 1 

u{dy) = 2b(y)dy + S (dy) 

a 

with b(y) = —A sgn(y) , y G E as in (6.1) above, and with S (-) the Dirac mass at the origin. 
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We shall deal with the equation (6.2) using a direct methodology that removes the parts of finite 
variation, that is, both the drift and the local time, and "reduces" (6.2) to a stochastic differential 
equation in natural scale 

Z(.) = p(y )+ [ s(Z(t))dW(t) (6.4) 
Jo 

for appropriate functions p(-) and s(-) . This approach was pioneered for the skew Brownian 

motion itself (i.e., with A = 0) by HARRISON & SHEPP HARRISON, J.M. & Shepp, L.A. 

(1981) (1981) , and for more general equations of the form (6.3) for suitable measurable functions 
r (•) and measures v on B(R) , by Le Gall Le Gall, J.F. (1983) (1983) Le Gall, J.F. (1984) 
(1984) and Engelbert & Schmidt Engelbert, H.J. & Schmidt, W. (1984) (1984). The 

results in these works do not seem to cover the equation (6.2), but those in BASS & Chen R. 

BASS, & Chen, Z.Q. (2005) (2005) do; we have preferred to detail a direct construction which 

is, in our opinion at least, quite simpler. 

In this spirit, let us introduce the scale function 

P(y) = (e 2X v - I) , y>0; p(0) = 0; p(y) = (l - e" 3 *) , y < . 

This has left-continuous derivative 

p\y) = (l-a)e 2Xy l { o,oo)(y) + ae- 2X yi { ^ ](y), y G R 

which is bounded away from zero, and second derivative measure 

p"(dy) = -2 b(y) p'(y) dy + (l - 2 a) S (dy) . 

Likewise, we introduce the inverse 

1 f 2\z \ , , -1 , / 2\z \ 

q(z) = — log [l + — J, ,>0; ,(0) = 0; qW = _log(l-— j, ,<0 

of the function p(-) , as well as its left-continuous derivative and its second-derivative measure 

q'O) = (1 - a + 2\z)~ 1 l (0jOo) (z) + (a - 2 A l ( _ OOj0 ](^) , 

q"(dz) = 2b(z) (q'(z)) 2 dz+ 2 "~\ 6 (dz) . 

a {1 — a) 

Analysis: Assume that a solution to (6.2) has been constructed; in particular, the process Y(-) 
is then a continuous semimartingale for which (3.14) holds a.s. We look at the process Z(-) : = 
p (¥(■)) and apply the Ito-Tanaka rule 

dZ(t) = p / (r(t))dF(t) + d / L Y (t,y)p"(dy) 

of (2.9), to obtain 

dZ(t) = p'(Y(t)) \b(Y(t))dt + dW(t)+ 2a ~ 1 dL Y (t) 
L a 
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-d / 2b(y)p'(y)L Y (t,y)dy + (l-2a)dL Y (t) = p' (q(Y (t))) dW (t) . 
Jr 

We have used here the occupation-time-density formula (2.8), and the property p'(0) = a . Now 
the piecewise-linear function 

s(z) := p'(q(z)) = = (l-a+2Xz)l {0<oo) (z)+(a-2Xz)l^ oo<Q] (z), zeR (6.5) 

9 \ z ) 

is bounded away from the origin, so the process Z(-) is the pathwise unique, strong solution of the 
stochastic differential equation (6.4) for this new dispersion function (Nakao Nakao, S. (1972) 
(1972)); and because Y(-) and Z(-) are bijections of each other, we have again the filtration 
identities 

3 z (t) = $ Y (t) = $ w (t) , < t < oo . (6.6) 

Synthesis: Consider the strong solution Z(-) of the stochastic differential equation (6.4) with 
the new dispersion function of (6.5), and define the process Y(-) := q(Z(-)) . Since dZ(t) = 
s(Z(t)) dW(t) and d(Z)(t) = s 2 (Z(t)) dt , this process satisfies almost surely 

roc rco d(Z)(t) 2 f°° 

l {z{ t)=o}dt = l{z(t)=o} 5 2(^( t )) - ( min («' i-"))" 2 J o Mz(t)=o}d(Z)(t) = 0, 
and we apply the ITO-TANAKA rule to obtain 

dY(t) = q'{Z{t))dZ(t)+d [ L z {t,z)q"{dz). 

Jr 

On the strength of the occupation-time-density formula and s(-) q'(-) = 1 , this gives 

dY(t) = dW(t)+d [ 2b(y) (q'(y)) 2 L Y (t,y)dy+ 2 * ~ 1 dL z {t) 
Jr « (1 - a) 

= dW(t) + b(Y(t)) {q'(Y(t))) 2 B 2 (Y(t))dt + ^~_ l a) dL z (t) 

= dW(t) + b(Y(t))dt + 2a ~ l dL Y (t), 

a 

that is, the equation (6.2) of the skew Brownian motion with bang-bang drift for the process Y(-) . 
We have used here the comparison of the local times at the origin for these two processes: 

L z {.) = (l-a)L Y (-). (6.7) 

This last identity (6.7) can be justified as follows: We start by noting 

LZ ^ = / l{o<z(t)< £ } d <^)(*) = lim T7 I l{o<y(t)<,( £ )} {p'(Y(t))) 2 dt 



l- a f fp'(Y(t))\ 
1 - a) -hm — / l{ <y(t)<q( £ )} I — : I dt . 



2 

Vlu' 2; ./o \ 1 - n ) 
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On the event {0 < Y(t) < q(e)} we have 

1 < p/(F(t)) <e 2 ^; and since lim f ^ 
1 — a e\.o \ e J 

we deduce the claimed identity of (6.7), namely 

L Z (-) = (1 - a) • Km ^ ^ l{o<y W < q ( £ ) } dt = (1 

• Taken together, the Analysis and Synthesis parts of this argument establish the following result. 

Theorem 6.1. The equation (3.18) admits a pathwise unique, strong solution for all values of its 
"skewness parameter" a G [0, 1] , and we have the filtration identity $ ¥ (t) = $ w (t) , < t < oo 
in (3.19). 

Remark 6.1. As shown towards the end of section 3 in Harrison & Shepp Harrison, J.M. 
& Shepp, L.A. (1981) (1981) , the stochastic equation (3.18) admits no solution for a [0, 1] . 
Consequently, when r] + £ ^ holds but the condition (1.5) fails because a — t]/(r) + £) ^ [0, 1] , 
the system of equations (1.2), (1.3) admits no solution. 

Remark 6.2. We compute in the next subsection the transition probabilities of the diffusion Y(-) . 
It follows from these computations, and in conjunction with the theory developed in Portenko 
Portenko, N.I. (1979) (1979) Portenko, N.I. (1979.b) (1979) Portenko, N.I. (1990) (1990) 
, that this process has the strong Markov and Feller properties. 

From the Remarks 3.1, 3.2, 6.1 and in conjunction with Theorem 4.1, we obtain now the 
following result. 

Proposition 6.1. The conditions of (1.5) are not just sufficient but also necessary for the well- 
posedness of the system of equations (1.2), (1.3). 

6.1 Joint Distribution of SBBBM and its Local Time 

Let us recall a construction of the skew Brownian motion (Ito & Mc Kean Ito, K. & Mc Kean, 
H.P, Jr. (1963) (1963), Ito, K. & McKean, H.P, Jr. (1974) (1974), Walsh Walsh, J.B. 
(1978) (1978)). We take a Brownian motion starting from y > 0, reflect it at the origin, and 
consider its excursions away from the origin. Then we change the sign of each excursion inde- 
pendently with probability 1 — a G (0, 1) . The resulting process is positive with probability a , 
negative with probability 1 — a . This implies a non-symmetric reflection principle around the 
origin. We shall see that even in the presence of the bang-bang drifts 

b(y) = -Asgnfo), y G R 

as in (6.1), this principle continues to hold for the skew Brownian motion. Thus, the joint distribu- 
tion of SBBBM and its local time are derived here. 

By GiRSANOV's theorem (e.g., Karatzas & Shreve Karatzas, I. & Shreve, S.E. 
(1991) (1991) , section 3.5), we consider the "reference probability measure" P* , under which the 



1 

1 — a ' 



-a)L Y {.). 
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process f Q b(Y(t)) dt + W(-) becomes standard Brownian motion. For every given t E [0, oo) 
the Radon-Nikodym derivative on # y (i) of the original measure with respect to the reference 
measure, is 



dP 



exp{ f b(Y(s))dW(s) + - [ b' 



(Y(s))ds 



exp{\(\y \-\Y(t)\+2L Y (t))-^t) 



we have used in this last equation the relationships (3.29), (3.9). 

Under the reference probability measure P* , the process Y(-) is skew Brownian motion which 
starting at y . As shown in Walsh Walsh, J.B. (1978) (1978) (see also Lang Lang, R. (1995) 
(1995)), the transition probability density function p^(t;y 0} ^) = F*(Y(t) E d£) /d£ for this 
process is given by 



P*(*; 3/o,0 



exp 



2t 



(2a 



sgn(0 



'2nt 



exp 



2t 



for (£, y) E M 2 , t > . Moreover, by the method of elastic Brownian motion (e.g., KARATZAS 
& Shreve Karatzas, I. & Shreve, S.E. (1984) (1984), Appuhamillage et al. Ap- 

PUHAMILLAGE, TH., VRUSHALI, B., THOMANN, E., WAYMIRE, E. & WOOD, B (2011) (2011)) 
the joint distribution of the skew Brownian motion and its symmetric local time is computed as 



P*(Y(t) E df , 2L Y (t) E db) = 
= {l + (2a-l)sgn(£)} 



\Z\+b+\y \ 
V2^ 



exp 



m+b+\ yo \y 

2t 



and 



F,(Y(t)Ed^2L Y (t) = 0) 



1 



2nt L 



exp 



(|yol - W 

2t 



exp 



Jd^dfe; 6>0 



2t 



d£ 



for £ E K . Note that we have 

1 + (2a - 1) sgn(0 = 2a if £ > , and 1 + (2a - 1) sgn(f) = 2(1 - a) if £ < . 

Thus, the non-symmetric reflection principle around the origin works intuitively even for the joint 
distribution. As expected, when there is no accumulation of local time at the origin, the skewness 
parameter a does not affect the transition probabilities. 

• We bring the above formulae from the reference measure P* back to the original measure P , 
by means of 

P(y ± (t) E A , Y T (t) = 0, 2L Y (t) EB) = (6.8) 
= exp {X\y \ - XH/2} • E p * [ exp {2L Y (t) - Y±(t)} ■ V ±{t)e ^ (t) = , 2 2^ )efl} 
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for (A, B) G B(R) x B(R) , t > . With (£, 6) G [0, oo) x (0, oo) , the joint density functions 



are 



¥(Y + (t) e d£,Y~(t) = 0, 2L Y (t) G d&) 



(6.9) 



2a ■ e 



-2A£ . £ + \yo\ 



exp 



+ + \y \ - At)' 

2t 



d^db. 



as well as 



P(Y~(t) G d£,Y + (t) = 0, 2L y (t) G d&) 
= 2(1 - a) -e 



(6.10) 



. i±A±N exp r _ (j+A+N^M: l de d6 . 



\/27t3 ~~~ x " I 2 t 

Whereas, when there is no accumulation of local time, we have 

F(Y ± (t) G d£,Y T (t) = 0, 2L Y (t) = 0) = 
1 



2ixt 



, (e-bol + At) 2 ) _ 2A5 f (£+M + At)< 
exp < > — e s • exp 1 



2t 



2t 



(6.11) 

})de, e>o. 



• The marginal density p(t;y ,^) d^ = P(Y(t) G d£) of Y(t) under the original probability 
measure P is obtained from 

P(Y(t) G d£) = P(Y(t) G d£ , 2L y (t) > 0) + P(Y(t) G d£ , 2L y (t) = 0) • l Uyo>0} , (6.12) 

where the second term takes care of the case when the local time is absent. If £ > and y > , 
the marginal density becomes 

(£ + y - At) 2 



p(t;y ,0 = (2«-l)e- 2A « 



1 



exp 



+ 



1 



exp 



'2nt * L 2t 

(e - i/o + At) 2 



(6.13) 



/27Tt * I- 2t 
whereas, if £ < and y < , this expression becomes 



+ (2a) / e ^^dw; 



'2ixt 



-fid 



p(t;y ,0 = (l-2a)e 2 ^-— = exp 



(-£ - |/o - At) 5 



+ 



exp 



27rt * L 2t 
(-£ + + At) 2 



/2vrt * I 2t 
If £ 2/o < , then this expression becomes 

p(t;y ,0 = {l + (2a-l)sgn(0}-- 



2fl-aV 



2nt 



(6.14) 



(u-\t) z 

e 2t dM . 



6-2/0 



-2A|C| 



exp 



/27Tt 

(KI + bol-At) 2 
2t 



(n-At) 2 

+ A / e 2t d-u 
[fl+lwl 



(6.15) 



where the term P(F(t) G d£ , 2L (t) = 0) in (6.12) is now equal to zero. 
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Remark 6.3. Letting t — Y oo in (6.13)-(6.15), we derive for the process Y(-) the stationary 
measure m(-) = / poo(£)d£ with the double-exponential probability density function 

Poo (0 := limp(t;y ,0 = a ■ (2A) e - 2A « • 1 {5>0} + (1 - a) ■ (2A) e 2 A « • 1 {5 < 0} . (6.16) 

It can be verified that (6.16) is the invariant distribution. Furthermore, it follows from the transition 
density (6.13)-(6.15) and the stationary distribution (6.16) that the following duality holds: 

/ 9(y)( [ m)p(t;y,m)poo(y)dy = [ /(£)( / g(y)p(t;^y)dy)p 00 (^, 
Jr v Jr 7 Jr K Jr 7 

for arbitrary bounded, measurable functions f,g,or equivalently 

[ g(y)E y [f(Y(t))]p 00 (y)dy = //(OE^(y(t))] Poo (Od£; *>0. (6.17) 

Here E y stands for the expectation under the measure F y induced by Y(-) which starts from 
l/Gl. Thus, under the probability measure 

Poo(-) := f F y (-)m(dy), 
Jr 

the process Y(-) is stationary; and moreover, given a fixed time T G (0, oo) , the time reversal 

Y(t) := Y(T — t) , 0<t<T (6.18) 

satisfies 

E F °°[f (Y(t ))---f n (Y(t n ))] = E r °°[f (Y{t n ))---f n {Y(t ))] (6.19) 

for every integer n 6 N , collection of time points = to < t\ < • • • < t n = T, and bounded, 
measurable functions / , . . . , f n . 

Remark 6.4. The infinitesimal generator of the process Y(-) may be defined formally by 

[£/](£) := ^r(e)-Asgn(0/ / (e) + 2(2«-l)/'(0^o(0, (6.20) 

for / e V := Cq°(IR) , where 5q(-) is the "Dirac delta function" at the origin. Here we use the 
parametrization for the symmetric local time L Y (•) . Let us denote formally the symmetric version 
of the density of m by 

PootO = (2a) A e~ A « ■ 1 {€>0} + 2(1 - a) A e A « ■ 1 {?<0} + A ■ 1 {5 = 0} . 

Then by direct calculation 

/ /(0 [£</](£) m(d£) = / g(0 [Cf}(0 m(dO , / [Cf] m(d£) = ; f,g£V. (6.21) 

</R JR </R 

Applying Theorem 2.3 of Fukushima & Stroock Fukushima, M. & Stroock, D. (1986) 
(1986) , we arrive at the same conclusion (6.19). 
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Remark 6.5. Let us define the time reversal of Y(-) as in (6.18). Following Pardoux Pardoux, 
E. (1986) (1986) and Petit Petit, F. (1997) (1997) , we may show that the time reversal is a 
solution of the stochastic equation 



Y(t) = Y(0) + W\t) + 2(1- 2a) L Y "(t) + 

+ jf sgn(y(s)) + ^ log p (T - s; y , Y(s)^j ds 



(6.22) 



for < t < T , where W^(-) is a standard Brownian motion with respect to the backwards 
filtration F y (-) generated by the time-reversed process Y(-) of (6.18), and 



L Y (t) := L Y (T) — L Y (T — t 



Y I 



0<t<T. 



(6.23) 



In the special case y — = Y(T) , the logarithmic derivative of the transition probability density 
function is 



where 

:= exp 



2t 



£ 2 (t,0 := Ae" 2A l«l / exp 



lei 



(u - At) : 
2t 



Thus, the time reversal is a skew Brownian bridge with bang-bang drift 
?{■) = Y(0) + W\-) + 2(l-2a)L Y (-)- 



\sgn(Y(t)) + 



Y(t) ( d 



(r-t,y(t)) 



This result suggests that the time reversal of SBBBM in general looks like a skew Brownian bridge 
drifted towards the target point Y(T) = y . 



7 Applications of the Skew Representations (3.30)-(3.31) 

With the skew representations in section 3.4 and the joint distribution of (Y (■) , L Y (■)) in section 
6.1 it is now straightforward to compute the transition density of the system (1.2)-(1.3) as well as 
its time reversal. 



7.1 Transition Density 

Let us discuss only some special cases, since the other cases are quite similar. For example, in the 
degenerate case with a = , thus p = 1 , 7 = 1 and with x% > %2 , (3.30)-(3.31) become 

X 1 (t) = x 1 + gt+ (Y + (t) - y+) - (3 L Y (t) , X 2 (t)=x 2 + gt+ (Y~ (t) - y ) - (3 L Y (t) , 
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for < t < oo , where y = x\ — x 2 > , and hence the transition density of (Xi(-), X 2 (-)) is 

HMt) e d6 , x 2 (t) e d£ 2 ) = (2a) • - ■ e -^-&) . exp { - ^^ foidfe , 



, 2 + /3\, f 2 ~P\ 2 

where Ci := £i - — - — £2 + ^1+ — 5 — p^ + ^-fl't. 



if > , £1 > £ 2 and £ 2 < x 2 + gt . Similarly, by (skew) symmetry: 

P(X 1 (t)Gd6, X 2 (t)Gd6) = 2(1 -a) • I • e- 2 ^^> ■ exp { - { ' 2 ~ A ^ }dfrd£ 

where c 2 := £ 2 - — 5 — Ki + % + — 5 — )x 2 + -z gt. 



2 



if /3 > 0, £ 2 > 5 and £1 < x 2 + gt. If f3 > , £1 > £ 2 = x 2 + g t , then the local time L Y '(•) 
is absent, and the transition density is easily obtained from (6.11): 

F(X 1 (t) edfi, X 2 (t) = x 2 + gt) = 

a - xi + x 2 + At) 2 1 _ 2Xa r (a + x\ - x 2 + At) 



exp < >— e -exp 



2irt V I 2t J I 2* 



d£i 

« = ?i-^2-9< 



• For another extreme example, in the degenerate case with a = 1 , thus p = 0, 7 = — 1 and 
with £1 > x 2 , (3.30)-(3.31) become 

X 1 (t) = x 1 -ht-{Y-{t)-yM2-0)L Y {t) , X 2 (t) = x 2 -ht-(Y + (t)-y+) + (2-f3) L Y (t) , 
for < t < 00 . If (3 < 2 , £1 > £ 2 and £1 > Xi - h t , then 

P(*i(f)edfc, X 2 (t)Gd6) = ( 2 ^)-^-e- 2A(€l - 6) - 7 ^exp{- (C3 ~ t At)2 }d£ 1 d6, 
If /3 < 2 , £ 2 > £1 and £ 2 > x 1 - ht , then 

F{X 1 (t)ed^ 1 , X 2 (t)Gd6) = 2(1 -a) —p --^=exp{ - U ^ > }dfrd&, 

If /3 < 2, £1 = Xi — ht > £ 2 , then the local time L y (-) does not accumulate, that is, the 
transitions density is obtained from (6.11) : 

P(Xi(t) = x x -/it, X 2 (t) G d£ 2 ) = 
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x 1 + x 2 + \t) 2 } _ 2Aq r (a + x 1 -x 2 + Xt) 



exp < > — e exp 



a = xi—fe — ht 



• For the isotropic variance case with p = a = 1 / \J2 ,7 = and xi > x 2 the difference and 
the sum of Xi(-) and X 2 (-) are 

X 1 (-) - X 2 {-) = ¥{■) , X x (-) + X 2 (-) = Xl + x 2 + vt + 2(2« - 1)L Y (-) + Q(-) , 

where (Y(-),L Y (-)) and Q(-) are independent. Thus the joint distribution of (Xi(-), X 2 (-)) are 
obtained by integrating out the local time. 

If a G (0, 1) \ {1 / 2} , then the above transition densities are discontinuous on the diagonal 
line due to the skewness. If a = 1 / 2 and (3 = 1, then these formulae are the same as those of 
the degenerate system studied in Fernholz et al. Fernholz, E.R., Ichiba, T., Karatzas, 
I. & Prokaj, V. (2012) (2012). The transition densities for all the other cases as well as the joint 
distribution of (Xx(-) , X 2 (-) , L Xl ~ X2 (-)) are computable from the skew representations (3.30)- 
(3.31) and the joint distribution (6.9)-(6. 11) in a similar manner. 

7.2 Time Reversal 

We consider now the time-reversal 

Xi(t) := Xi(T — t) , Xi{t) := Xi(T-t)-Xi(T), 0<t<T,i = l,2 (7.1) 

of the solution to the system (1.2)-( 1.3) with the backwards filtration F = {d(t)}o<t<T generated 
by the random variable Y(T) and by the time-reversal 

(W(t) := W(T-t)-W(T), Q(t) := Q(T - t) - Q(T) ) , 0<t<T 

of the planar Brownian motion : 

$(t) := *(Y(T)) \/$ im \t), $& W \t) := <t(Q(0),W(0) ;0 < < t) , 0<t<T. 

With some extra work in addition to the discussion of Remark 6.5 we may show that Y(-) is a 
diffusion (6.22) driven by the F -Brownian motion W\-) (c.f. Pardoux Pardoux, E. (1986) 
(1986) and section 3 of Petit Petit, F. (1997) (1997)). Combining the skew representations 
(3.30)-(3.3 1) with the time-reversals (6.22)-(6.23), we derive the time-reversed skew representation 
form 

Mt) = -^t + p 2 (r+(t)-F+(0))-a 2 (F-(t)-F-(0))-(l-/3- 7 )2 y (t)+p ( Tg(t) 

= -lit+ [ {p 2 l {H +a 2 l {H )dY(s)-(l-(3-2 1 )L Y (t)+paQ(t) (7.2) 
for < t < T , and 

X 2 (t) = -pt-a 2 (F + (t)-r + (0))+p 2 (F-(t)-F-(0))-(l-/3- 7 )2 y (t)+P^g(t) 

= —fit — j {p 2 l {? +a 2 l {? y\Y{s)-(l-P-2 1 )L Y (t)+poQ(t). (7.3) 
^0 
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Remark 7 .1. By analogy with (7.1) we denote the time-reversal of ranks by Ri(t) := Ri(T — t) 
for < t < T, i — 1, 2 . Applying the Tanaka formula to (7.2)-(7.3), we may derive the 
time-reversed dynamics of (-Ri(-), R2(-)) • 

Remark 7.2. As we saw in Remarks 6.3-6.4, the process Y(-) is strictly time-reversible when 
started at its invariant distribution (6.16). Under this invariant distribution, the dynamics of the 
time-reversal of (Ri(-), Rz(')) can De derived through the skew representations of (3.30)-(3.31). 
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